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Strongly correlated Fermi systems are among the most intriguing, best experimentally studied 
and fundamental systems in physics. These are, however, in defiance of theoretical understanding. 
The ideas based on the concepts like Kondo lattice and involving quantum and thermal fluctuations 
at a quantum critical point have been used to explain the unusual physics. Alas, being suggested to 
describe one property, these approaches fail to explain the others. This means a real crisis in theory 
suggesting that there is a hidden fundamental law of nature, which remains to be recognized. A 
theory of fermion condensation quantum phase transition, preserving the extended quasiparticles 
paradigm and intimately related to unlimited growth of the effective mass as a function of tem- 
perature, magnetic field etc, is capable to resolve the problem. We discuss the construction of the 
theory and show that it delivers theoretical explanations of vast majority of experimental results in 
strongly correlated systems such as heavy-fermion metals and quasi-two-dimensional Fermi systems. 
Our analysis is placed in the context of recent salient experimental results. Our calculations of 
the non-Fermi liquid behavior, of the scales and thermodynamic and transport properties are 
in good agreement with the heat capacity, magnetization, longitudinal magnetoresistance and 
magnetic entropy obtained in remarkable measurements on the heavy fermion metal YbRli2Si2. 
Using two-dimensional 3 He as an example, we demonstrate that the main universal features of its 
experimental temperature T - density x phase diagram resemble those of the heavy-fermion metals. 
We propose a simple expression for the effective mass, describing all diverse experimental facts 
on the 3 He in unified manner and demonstrating that the universal behavior of the effective mass 
coincides with that observed in heavy fermion metals. 
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I. INTRODUCTION 

Strongly correlated Fermi systems represented by 
heavy fermion (HF) metals, high-temperature supercon- 
ductors and quasi- two-dimensional 3 He are among the 
most intriguing, best experimentally studied and funda- 
mental systems in physics This is also a field never 
far from applications in synthesis of novel materials for 
cryogenics, rare earth magnets and applied superconduc- 
tivity. Their behavior is so unusual that the traditional 
Landau quasiparticles paradigm does not apply to it. 
The paradigm states that the properties is determined 
by quasiparticles whose dispersion is characterized by the 
effective mass M* which is independent of temperature 
T, the number density x, magnetic field B and other ex- 
ternal parameters. The above systems are, however, in 
defiance of theoretical understanding. The ideas based 
on the concepts (like Kondo lattice, see e.g. Ref. @) in- 
volving quantum and thermal fluctuations at a quantum 
critical point (QCP) have been used to explain the un- 
usual physics of these systems known as non- Fermi liquid 
(NFL) behavior [JQ. Alas, being suggested to describe 
one property, these approaches fail to explain the others. 
This means a real crisis in theory suggesting that there 



'Electronic address: vrshag@thd.pnpi.spb.ru 



is a hidden fundamental law of nature, which remains to 
be recognized [f|. It is widely believed that utterly new 
concepts are required to describe the underlying physics. 
There is a fundamental question: how many concepts do 
we need to describe the above physical mechanisms? This 
cannot be answered on purely experimental or theoretical 
grounds. Rather, we have to use both of them. 

Usual arguments that quasiparticles in strongly corre- 
lated Fermi liquids "get heavy and die" at a quantum 
critical point commonly employ the well-known formula 
basing on assumptions that the z-factor (the quasipar- 
ticle weight in the single-particle state) vanishes at the 
points of second-order phase transitions [6|. However, it 
has been shown that this scenario is problematic (7j. A 
concept of fermion condensation quantum phase transi- 
tion (FCQPT) preserving quasiparticles and intimately 
related to the unlimited growth of M* , had been sug- 
gested |§-HH- Studies show that it is capable to deliver 
an adequate theoretical explanation of vast maj ority of 
experimental results in different HF metals [13 - il4| . In 
contrast to the Landau paradigm based on the assump- 
tion that M* is a constant, in FCQPT approach M* 
strongly depends on T, x, B etc. Therefore, in accord 
with numerous experimental facts the extended quasi- 
particles paradigm is to be introduced. The main point 
here is that the well-defined quasiparticles determine as 
before the thermodynamic and transport properties of 
strongly correlated Fermi-systems, while M* becomes a 
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function of T, x, B, and the dependence of the effective 
mass on T, x, B gives rise to the non-Fermi liquid (NFL) 
behavior pj3, Mu^- 

In this review report we discuss the construction of a 
theory, based on the above FCQPT approach and its ap- 
plication to the analysis of wide variety of experimentally 
observed phenomena in microscopically different strongly 
correlated Fermi systems like heavy-fermion metals and 
quasi-two-dimensional 3 He. We analyze the NFL behav- 
ior of strongly correlated Fermi systems and show that 
this is generated by the dependence of the effective mass 
on temperature, number density and magnetic field at 
FCQPT. We demonstrate that the NFL behavior ob- 
served in the transport and thermodynamic properties 
of HF metals can be described in terms of the scaling be- 
havior of the normalized effective mass. This allows us to 
construct the scaled thermodynamic and transport prop- 
erties extracted from experimental facts in wide range of 
the variation of scaled variable. We show that "pecu- 
liar points" of the normalized effective mass give rise to 
the energy scales observed in the thermodynamic and 
transport properties of HF metals. Our calculations of 
the thermodynamic and transport properties are in good 
agreement with the heat capacity, magnetization, lon- 
gitudinal magnetoresistance and magnetic entropy ob- 
tained in remarkable measurements on the heavy fermion 
metal YbRh 2 Si 2 fTHHI. 



II. 



FERMION CONDENSATION QUANTUM 
PHASE TRANSITION 



We start with visualizing the main properties of FC- 
QPT. To this end, consider the density functional theory 
for superconductors (SCDFT) [HJ]. SCDFT states that 
at fixed temperature T the thermodynamic potential $ is 
a universal functional of the number density n(r) and the 
anomalous density (or the order parameter) «(r,ri) and 
provides a variational principle to determine the densi- 
ties [22[ . At the superconducting transition temperature 
T c a superconducting state undergoes the second order 
phase transition. Our goal now is to construct a quantum 
phase transition which evolves from the superconducting 
one. In that case, the superconducting state takes place 
at T = while at finite temperatures there is a normal 
state. This means that in this state the anomalous den- 
sity is finite while the superconducting gap vanishes. For 
the sake of simplicity, we consider a homogeneous Fermi 
(electron) system. 

Let us assume that the coupling constant A of the pair- 
ing interaction vanishes, A — > 0, making vanish the su- 
perconducting gap at any finite temperature. In that 
case, T c — > and the superconducting state takes place 
at T = while at finite temperatures there is a normal 
state. This means that at T = the anomalous density 
is finite while the superconducting gap is infinitely small 
[Tol . Il2l . [23j j . For the sake of simplicity, we consider a ho- 
mogeneous Fermi (electron) system [1 21 ] . Then, the ther- 



modynamic potential reduces to the ground state energy 
E which turns out to be a functional of the occupation 



number n(p) since K — \/n(l — n) [I6|, |17|, |22J, |24|, [2 
Upon minimizing E with respect to n(p), we obtain 
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where /i is the chemical potential. It is seen from Eq. ([T]) 
that instead of the Fermi step, we have < n(p) < I in 
certain range of momenta pi < p < pf with k is finite in 
this range. Thus, the step-like Fermi filling inevitably un- 
dergoes restructuring and formes the fermion condensate 
(FC) as soon as Eq. ([T]) possesses not-trivial solutions at 
some point x — x c when pi = pj = pp [H, [T^, EH- Here 
Pf is the Fermi momentum and x = p F /3ir 2 . 

At any small but finite temperature the anomalous 
density k (or the order parameter) decays and this state 
undergoes the first order phase transition and converts 
into a normal state characterized by the thermodynamic 
potential $ - At T — > 0, the entropy S = -<9<J> /<9T of 
the normal state is given by the well-known relation (2(| 



S = -2 / [n(p) ln(n(p)) + (I - n(p) ln(l - n(p))] 



dp 

(27r)3' 

which follows from combinatorial reasoning. Since the 
entropy of the superconducting ground state is zero, it 
follows from Eq. ([2]) that the entropy is discontinuous at 
the phase transition point, with its discontinuity AS = 
Sq. The latent heat q of transition from the asymmetrical 
to the symmetrical phase is q = T c Sq — since T c = 
0. Because of the stability condition at the point of the 
first order phase transition, we have $o["(p)] = $[ K (p)]- 
Obviously the condition is satisfied since q = 0. 

At T = 0, a quantum phase transition is driven by a 
nonthermal control parameter, e.g. the number density 
x. To clarify the role of x, consider the effective mass M* 
which is related to the bare electron mass m by the well- 
known Landau equation [2(| which is also valid when M* 
strongly depends on B, T or x Hi: 
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Here we omit the spin indices for simplicity, n(p, T) is 
quasiparticle occupation number, and F is the Landau 
amplitude. At T = 0, Eq. © reads H^Hl 



M* 
rn 



1 



1 - N F 1 (x)/3 ' 



(4) 



Here Nq is the density of states of free electron gas and 
F (x) is the p-wave component of Landau interaction 
amplitude F. When at some quantum critical point 
(QCP) x = x c , F 1 (x) achieves certain threshold value, 
the denominator in Eq. ^ tends to zero so that the effec- 
tive mass diverges at T — [13, HI] . It follows from Eq. 
(j4|) that beyond the QCP x = x c , the effective mass be- 
comes negative. To avoid unstable and physically mean- 
ingless state with a negative effective mass, the system 
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must undergo a quantum phase transition at QCP x = x c 
denned by Eq. Q and which is FCQPT % i, El El- 
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FIG. 1: Schematic phase diagram of the system driven to 
the FC state. The number density x is taken as the control 
parameter and depicted as x/x c . The quantum critical point 
(QCP), x/x c = 1, of FCQPT is shown by the arrow. At 
x/x c < 1 and sufficiently low temperatures, the system is in 
the Landau Fermi liquid (LFL) state as shown by the shadow 
area. At T — and beyond QCP, x/x c > 1, the system is at 
the quantum critical line depicted by the dash line and shown 
by the vertical arrow. The critical line is characterized by the 
FC state with finite superconducting order parameter k. At 
any finite temperature, the order parameter k is destroyed 
and the entropy becomes discontinuous at T c = 0, the system 
undergoes the first order phase transition and exhibits the 
NFL behavior at T > 0. 

Schematic phase diagram of the system which is driven 
to FC by variation of x is reported in Fig. [T] Upon ap- 
proaching the critical density x c the system remains in 
LFL region at sufficiently low temperatures [H El| , that 
is shown by the shadow area. At QCP x c shown by the 
arrow in Fig. [TJ the system demonstrates the NFL be- 
havior down to the lowest temperatures. Beyond QCP 
at finite temperatures the behavior is remaining the NFL 
one and is determined by the temperature-independent 
entropy So (HI- In that case at T — ¥ 0, the system is ap- 
proaching a quantum critical line (shown by the vertical 
arrow and the dashed line in Fig. QJ rather than a quan- 
tum critical point. Upon reaching the quantum critical 
line from the above at T — > the system undergoes the 
first order quantum phase transition, which is FCQPT 
taking place at T c = 0. 

At T > the NFL state above the critical line, see Fig. 
[TJ is strongly degenerated, therefore it is captured by the 
other states such as superconducting (for example, by 
the superconducting state in CeCoIns I23L |29| . 30] ) or by 
AF state (e.g. AF one in YbRh 2 Si 2 [SHOl) lifting 
the degeneracy. The application of magnetic field B > 
B c o restores the LFL behavior, where B c q is a critical 
magnetic field, such that at B > B c q the system is driven 



towards its Landau Fermi liquid (LFL) regime [U, |3(| 
l3lj ]. In some cases, for example in HF metal CcPai2Si2, 
fU = 0, see e.g. [32|, while in YbRh 2 Si 2 , B c0 ~ 0.06 T 
[2lll3~l|. In our simple model B c q is taken as a parameter. 



III. SCALING BEHAVIOR OF THE EFFECTIVE 
MASS 

Schematic phase diagram of the HF metal YbRh2Si2 
is reported in Fig. [5] Magnetic field B is taken as the 
control parameter. The FC state and the region lying 
at x/x c > 1, see Fig. [JJ can be captured by the super- 
conducting, ferromagnetic, antiferromagnetic (AF) etc. 
states lifting the degeneracy [H EI| ■ Since we consider 
the HF metal YbRh 2 Si 2 the AF state takes place HHU 
as shown in Fig. [5J As seen from Fig. [2J at elevated 
temperatures and fixed magnetic field the NFL regime 
occurs, while rising B again drives the system from NFL 
region to LFL one. Below we consider the transition re- 
gion when at rising B the system moves from NFL regime 
to LFL one along the dash-dot horizontal arrow, and at 
elevated T it moves from LFL regime to NFL one along 
the solid vertical arrow. 




Control parameter, magnetic field B 



FIG. 2: Schematic phase diagram of the HF metal YbRli2Si2. 
AF denotes antiferromagnetic state. At T = 0, B c o is mag- 
netic field at which the effective mass diverges and the AF 
state vanishes. At B > B c o the system is in its paramagnetic 
state. The vertical arrow shows the transition from the LFL 
regime to the NFL one at fixed B along T with M* depending 
on T. The dash-dot horizontal arrow illustrates the system 
moving from NFL regime to LFL one along B at fixed T. 
The inset shows a schematic plot of the scaling behavior of 
the normalized effective mass versus the normalized temper- 
ature. Transition regime, where Mjy reaches its maximum 
value at T — Tm, is shown by the hatched area both 
in the main panel and in the inset. The arrows mark the 
position of inflection point in Mjy and the transition region. 
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To explore a scaling behavior of M* , we write the 
quasiparticle distribution function as ni(p) = n(p,T) — 
n(p), with n(p) is the step function, and Eq. © then 
becomes 
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(5) 

At QCP the effective mass M* diverges and Eq. © 
becomes homogeneous determining M* as a function of 
temperature 



M*(T) oc T 



-2/3 



(6) 



while the system exhibits the NFL behavior [l2|, HH • If 
the system is located before QCP, M* is finite, at low 
temperatures the system demonstrates the LFL behav- 
ior that is M*(T) ~ M*+aiT 2 1 with a\ is a constant, see 
the inset to Fig. [2] Obviously, the LFL behavior takes 
place when the second term on the right hand side of Eq. 
([5]) is small in comparison with the first one. Then, at ris- 
ing temperatures the system enters the transition regime: 
M* grows, reaching its maximum at T = Tm, with 
subsequent diminishing. Near temperatures T > Tm the 
last "traces" of LFL regime disappear, the second term 
starts to dominate, and again Eq. © becomes homoge- 
neous, and the NFL behavior restores, manifesting itself 
in decreasing M* as T~ 2 / 3 , see Eq. ([6]). When the sys- 
tem is near QCP, it turns out that the solution of Eq. 
© M*(T) can be well approximate d by a simple uni- 
versal interpolating function [TJ|, [HI l30| . The interpo- 
lation occurs between the LFL (M* ~ M* + aiT 2 ) and 
NFL [M* oc T -2 / 3 ) regimes thus describing the above 
crossover [l2|, [l5[ . Introducing the dimensionless variable 
y = Tn = T /Tm, we obtain the desired expression 



M* N {y)^c 



1 + ay 2 



1 



ay 



8/3 ■ 



(7) 



Here = M* /M M is the normalized effective mass, 
Co = (1 + C2)/(l + Ci), ci and C2 are fitting parameters, 
parameterizing the Landau amplitude. 

The inset to Fig. [5] demonstrates the scaling behavior 
of the normalized effective mass Mjy = M*/M M ver- 
sus normalized temperature Tjv = T /Tm, where M M is 
the maximum value that M* reaches at T = Tm- At 
T < T M the LFL regime takes place. At T > T M 
the T -2 / 3 regime takes place. This is marked as NFL 
one since the effective mass depends strongly on tem- 
perature. The temperature region T ~ Tm signifies the 
transition between the LFL regime with almost constant 
effective mass and NFL behavior, given by T" 2 / 3 de- 
pendence. Thus temperatures T ~ Tm can be regarded 
as the transition region between LFL and NFL regimes. 
The transition temperatures are not really a phase tran- 
sition. These necessarily are broad, very much depending 
on the criteria for determination of the point of such a 
transition, as it is seen from the inset to Fig. [2] As 
usually, the transition temperature is extracted from the 



temperature dependence of charge transport, for exam- 
ple, from the resistivity p{T) = p Q + AT 2 with p is the 
residual resistivity and A is the LFL coefficient. The 
crossover takes place at temperatures where the resis- 
tance starts to deviate from the LFL T 2 behavior. Ob- 
viously, the measure of the deviation from the LFL T 2 
behavior cannot be defined unambiguously. Therefore, 
different measures produce different results. 

It is possible to transport Eq. (l5l) to the case of the 
application of magnetic fields [lU, [lfl H(| • The applica- 
tion of magnetic field restores the LFL behavior so that 



M M depends on B as 



M 



M oc {B — B, 
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while 



Tm oc ps(B - Bco), 



(8) 



(9) 



where ps is the Bohr magneton [I2|, |15|, |30(. Employ- 
ing Eqs. © and © to calculate M* M and T M , we 
conclude that Eq. ([7]) is valid to describe the normal- 
ized effective mass in external fixed magnetic fields with 
y = T/{B - B c0 ). On the other hand, Eq. © is valid 
when the applied magnetic field becomes a variable, while 
temperature is fixed T = Tt. In that case, as seen from 
Eqs. ([6]), (7} and©, it is convenient to rewrite both the 
variable as y = (B — B c o)/Tf, and Eq. © as 



Pb{B m - B c0 ) oc Tf. 



(10) 



It follows from Eq. © that in contrast to the Landau 
paradigm of quasiparticlcs the effective mass strongly de- 
pends on T and B. As we will see it is this dependence 
that forms the NFL behavior. It follows also from Eq. © 
that a scaling behavior of M* near QCP is determined 
by the absence of appropriate external physical scales to 
measure the effective mass and temperature. At fixed 
magnetic fields, the characteristic scales of temperature 
and of the function M* (T, B) are defined by both Tm and 
M M respectively. At fixed temperatures, the character- 
istic scales are (£>m — B c q) and M M . It follows from Eqs. 
© and © that at fixed magnetic fields, Tm —> 0, and 
M M — > 00, and the width of the transition region shrinks 
to zero as B —¥ B c q when these are measured in the exter- 
nal scales. In the same way, it follows from Eqs. © and 
(fTU)) that at fixed temperatures, (Bm — B c o) — > 0, and 
M M —7> 00, and the width of the transition region shrinks 
to zero as Tf — > 0. Thus, the application of the external 
scales obscure the scaling behavior of the effective mass 
and of the thermodynamic and transport properties. 

A few remarks are in order here. As we shall see, 
magnetic field dependencies of the effective mass or of 
other observable like the longitudinal magnetoresistance 
do not have "peculiar points" like maximum. The nor- 
malization are to be performed in the other points like 
the inflection point at T — Ti n f (or at B = Bi n f) shown 
in the inset to Fig. H]by the arrow. Such a normalization 
is possible since it is established on the internal scales, 
T in f oc T M oc [B - B c0 ). 
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IV. NFL BEHAVIOR OF THE HF METAL 
YbRh 2 Si 2 

In what follows, we compute the effective mass and 
employ Eq. ([7]) for estimations of considered values. To 
compute the effective mass M*(T,B), we solve Eq. §5§ 
with special form of Landau interaction amplitude, see 
Refs. [ill, [l5| for details. Choice of the amplitude is dic- 
tated by the fact that the system has to be at QCP, which 
means that first two p-derivatives of the single-particle 
spectrum e(p) should equal zero. Since first derivative 
is proportional to the reciprocal quasiparticle effective 
mass 1/M*, its zero just signifies QCP of FCQPT. Zero 
of the second derivative means that the spectrum e(p) 
has an inflection point at pf rather than a maximum. 
Thus, the lowest term of the Taylor expansion of e(p) is 
proportional to {p — Pf) 3 [HI- After solution of Eq. (JSJ, 
the obtained spectrum had been used to calculate the en- 
tropy S(B,T), which, in turn, had been recalculated to 
the effective mass M* (T, B) by virtue of well-known LFL 
relation M*(T, B) = S(T, B)/T. Our calculations of the 
normalized entropy as a function of the normalized mag- 
netic field B/Bi n f — y and as a function of the normal- 
ized temperature y = T/Ti n f are reported in Fig. [3] Here 
Ti n f and Bi n f are the corresponding inflection points in 
function S. We normalize the entropy by its value at 
the inflection point Sn(v) = S(y)/S(l). As seen form 
Fig. [31 our calculations corroborate the scaling behavior 
of the normalized entropy, that is the curves at different 
temperatures and magnetic fields merge into single one 
in terms of the variable y. The inflection point Ti n f in 
S(T) makes M* (T, B) have its maximum as a function of 
T, while M*(T,B) versus B has no maximum. We note 
that our calculations of the entropy confirm the validity 
of Eq. ([7]) and the scaling behavior of the normalized 
effective mass. 



A. Heat capacity 

Exciting measurements of C/T oc M* on samples of 
the new generation of YbRli2Si2 in different magnetic 
fields B up to 1.5 T [l9[ allow us to identify the scaling 
behavior of the effective mass M* and observe the dif- 
ferent regimes of M* behavior such as the LFL regime, 
transition region from LFL to NFL regimes, and the NFL 
regime itself. A maximum structure in C/T oc M M at 
Tm appears under the application of magnetic field B 
and Tm shifts to higher T as B is increased. The value 
of C/T = 70 is saturated towards lower temperatures 
decreasing at elevated magnetic field, where 70 is the 
Sommerfeld coefficient 19]. 

The transition region corresponds to the temperatures 
where the vertical arrow in the main panel of Fig. [2] 
crosses the hatched area. The width of the region, being 
proportional to Tm oc (B — B c q) shrinks, Tm moves to 
zero temperature and 70 oc M* increases as B — > B c q. 
These observations are in accord with the facts . 
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FIG. 3: The normalized entropy Sn(B / B in f) versus y — 
B/Bi n f an d the normalized entropy Sn(T /Ti n f) versus y = 
T/Ti n f calculated at fixed temperature and magnetic field, 
correspondingly, are represented by the solid lines and shown 
by the arrows. The inflection point is depicted by the dash- 
dot arrow. 



FIG. 4: The normalized effective mass extracted from 
the measurements of the specific heat C/T on YbRli2Si2 in 
magnetic fields B listed in the legend. Our calculations 
are depicted by the solid curve tracing the scaling behavior of 
M* N . 

To obtain the normalized effective mass M^, the max- 
imum structure in C /T was used to normalize C /T, and 
T was normalized by Tm- In Fig. H]M^ as a function of 
normalized temperature Tjv is shown by geometrical fig- 
ures, our calculations are shown by the solid line. Figure 
U reveals the scaling behavior of the normalized exper- 
imental curves - the scaled curves at different magnetic 
fields B merge into a single one in terms of the normal- 
ized variable y = T/Tm- As seen, the normalized mass 
extracted from the measurements is not a constant, 
as would be for LFL. The two regimes (the LFL regime 
and NFL one) separated by the transition region, as de- 
picted by the hatched area in the inset to Fig. [2J are 
clearly seen in Fig. U illuminating good agreement be- 



6 



tween the theory and the facts. It is worthy of note 
that the normalization procedure allows us to construct 
the scaled function C/T extracted from the facts in wide 
range variation of the normalized temperature. Indeed, it 
integrates measurements of C /T taken at the application 
of different magnetic fields into unique function which 
demonstrates the scaling behavior over three decades in 
normalized temperature as seen from Fig. 01 



B. Magnetization 

Consider now the magnetization M as a function of 
magnetic field B at fixed temperature T = Tf 



M(B,T) 



x(b,T)db, 



where the magnetic susceptibility X is given by 

j3M* (B, T) 



X(B,T) 



1 



F " 



(11) 



(12) 



Here, j3 is a constant and Fft is the Landau amplitude 
related to the exchange interaction (26[. In the case of 
strongly correlated systems F a > -0.9 [13, HI- There- 
fore, as seen from Eq. ([12"]) , due to the normalization the 
coefficients f3 and (1 + Ffi) drops out from the result, and 
XocM*. 



magnetization demonstrates the same behavior [2(j. We 
use Bk and M(Bk) to normalize B and M respectively. 
The normalized magnetization M(B)/M(Bk) extracted 
from facts [2(| depicted by the geometrical figures and 
calculated magnetization shown by the solid line are re- 
ported in Fig. [5J As seen, the scaled data at different Tf 
merge into a single one in terms of the normalized vari- 
able y = B/Tk- It is also seen, that these exhibit energy 
scales separated by kink at the normalized magnetic field 
Bn = B /Bk = 1 • The kink is a crossover point from the 
fast to slow growth of M at rising magnetic field. It is 
seen from Fig. [SJ that our calculations are in good agree- 
ment with the facts, and all the data exhibit the kink 
(shown by the arrow) at Bm ~ 1 taking place as soon as 
the system enters the transition region corresponding to 
the magnetic fields where the horizontal dash-dot arrow 
in the main panel of Fig. [2] crosses the hatched area. In- 
deed, as seen from Fig. at lower magnetic fields M is a 
linear function of B since M* is approximately indepen- 
dent of B. Then, it follows from Eqs. © and © that 
at elevated magnetic fields M* becomes a diminishing 
function of B and generates the kink in M[B) separat- 
ing the energy scales discovered in Refs. [Rj|2Q|]. Then, 
as seen from Eq. (fT0|) the magnetic field Bk at which the 
kink appears, Bk — Bm oc Tf, shifts to lower B as Tf is 
decreased. This observation is in accord with the facts 
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C. Longitudinal magnetoresistance 
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FIG. 5: The field dependencies of the normalized magneti- 
zation M collected at different temperatures shown at right 
bottom corner are extracted from measurements collected on 
YbRu2Si2 [20I ]. The kink (shown by the arrow) is clearly seen 
at the normalized field Bn = B/Bk — 1. The solid curve 
represents our calculations. 

One might suppose that Ftf can strongly depend on 
B. This is not the case, since the Kadowaki- Woods ratio 
is conserved [HHM1], A(B)/-yg(B) « A(B)/ X 2 (B) oc 
const, we have 70 oc M* oc x- Here A is the coefficient 
in the T 2 dependence of resistivity p. 

Our calculations show that the magnetization exhibits 
a kink at some magnetic field B = Bk ■ The experimental 
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0,01 0,1 1 10 

Normalized magnetic field 

FIG. 6: Magnetic field dependence of the normalized mag- 
netoresistance pM versus normalized magnetic field, pjv was 
extracted from LMR of YbRfi2Si2 at different temperatures 
18, 20] listed in the legend. The inflection point is shown by 
the arrow, and the solid line represents our calculations. 

Consider a longitudinal magnetoresistance (LMR) 
p(T, B) = po + A(T, B)T 2 as a function of B at fixed 
T = Tf. In that case, the classical contribution to 
LMR due to orbital motion of carriers induced by the 
Lorentz force is small, while the Kadowaki- Woods rela- 
tion [HI, HUSH, K = Aj-j 2 oc A/x 2 = const, allows us 
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to employ M* to construct the coefficient A [36], since 
7o oc x oc M*. As a result, p(T, B) - p oc (M*) 2 . Fig. EO 
reports the normalized magnetoresistance 

PN (y) = PM^L oc (M^(y)) 2 (13) 

Pin/ 

versus normalized magnetic field y = B/Bi n f at differ- 
ent temperatures, shown in the legend. Here pi n f and 
Bi n f are LMR and magnetic field respectively taken at 
the inflection point marked by the arrow in Fig. [6] Both 
theoretical (shown by the solid line) and experimental 
(marked by the geometrical figures) curves have been nor- 
malized by their inflection points, which also reveals the 
scaling behavior - the scaled curves at different temper- 
atures merge into single one as a function of the variable 
y and show the scaling behavior over three decades in 
the normalized magnetic field. The transition region at 
which LMR starts to decrease is shown in the inset to 
Fig. [5] by the hatched area. Obviously, as seen from Eq. 
(|10p . the width of the transition region being propor- 
tional to Bm — Bi n f oc Tf decreases as the temperature 
Tf is lowered. In the same way, the inflection point of 
LMR, generated by the inflection point of M* shown in 
the inset to Fig. [5] by the arrow, shifts to lower B as 
Tf is decreased. All these observations are in excellent 
agreement with the facts [13, |20j . 




T[K] 

FIG. 7: MR versus temperature T as a function of magnetic 
field B. The experimental data on MR were collected on 
CeCoIns at fixed magnetic field B [33] shown in the right 
bottom corner of the Figure. The solid lines represent our 
calculations. 

It is instructive to demonstrate that the same effective 
mass employed to calculate LMR shown in Fig. [5] gives 
good description of the magnetoresistance (MR) collected 
in measurements on CeCoIns. Figure [7] shows the calcu- 
lated MR versus temperature as a function of magnetic 
field B together with the experimental points from Ref. 
(37j . We note that both the classical contribution to MR 
due to orbital motion of carriers induced by the Lorentz 
force and po were omitted. As seen from Fig. our 
description of experiment is pretty good [38] . 



D. Magnetic entropy 

The evolution of the derivative of magnetic entropy 
dS(B,T)/dB as a function of magnetic field B at fixed 
temperature Tf is of great importance since it allows us 
to study the scaling behavior of the derivative of the ef- 
fective mass TdM*(B,T)/dB oc dS(B,T)/dB. While 
the scaling properties of the effective mass M* {B, T) can 
be analyzed via LMR, see Fig. [5] 




~i — ' — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — 1 — i — 1 — 
0123456789 10 

NORMALIZED MAGNETIC FIELD 



FIG. 8: Normalized magnetization difference divided by tem- 
perature increment (AM/AT)jv versus normalized magnetic 
field at fixed temperatures listed in the legend is extracted 
from the facts collected on YbRh2Si2 [2l[. Our calculation 
of the normalized derivative (dS/dB)N — (AM / AT)n versus 
normalized magnetic field is shown by the solid line. 

As seen from from Eqs. and (fTDj) . at y < 1 the 
derivative -dM N (y)/dy oc y with y = (B - B c0 )/(B in f - 
B c0 ) oc (B — B c0 )/Tf. We note that the effective mass as 
a function of B does not have the maximum. At elevated 
y the derivative —dM^(y)/dy possesses a maximum at 
the inflection point and then becomes a diminishing func- 
tion of y. Upon using the variable y = (B — B c o)/Tf, 
we conclude that at decreasing temperatures, the lead- 
ing edge of the function —dS/dB oc —TdM*/dB becomes 
steeper and its maximum at (B in f — B cQ ) oc Tf is higher. 
These observations are in quantitative agreement with 
striking measurements of the magnetization difference di- 
vided by temperature increment, — AM/ AT, as a func- 
tion of magnetic field at fixed temperatures Tf collected 
on YbRli2Si2 [H|. We note that according to the well- 
know thermodynamic equality dM/dT = dS/dB, and 
AM/ AT ~ dS/dB. To carry out a quantitative analy- 
sis of the scaling behavior of —dM*(B,T)/dB, we cal- 
culate the normalized entropy S shown in Fig. [3] as a 
function of B /Bi n f at fixed temperature Tf . Fig. [5] re- 
ports the normalized (dS/dB)^ as a function of the nor- 
malized magnetic field. The scaled function (dS/dB)^ 
is obtained by normalizing (~dS / dB) by its maximum 
taking place at Bm, and the field B is scaled by Bm- 
The measurements of — AM / AT are normalized in the 
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same way and depicted in Fig. [8] as (AM/ AT) n versus 
normalized field. It is seen from Fig. [5] that our calcula- 
tions shown by the solid line are in good agreement with 
the facts and the scaled functions (AAf/AT)jv extracted 
from the facts show the scaling behavior in wide range 
variation of the normalized magnetic field B/Bm- 

E. Energy scales and T — B phase diagram for 
YbRh 2 Si 2 



as well as these of LMR and the magnetic entropy are 
shifted to very low temperatures. Therefore due to ex- 
perimental difficulties these regimes cannot be often ob- 
served in experiments on HF metals. As it is seen from 
Figs. 21 [5l [U |8] and HI the normalization allows us to con- 
struct the unique scaled thermodynamic and transport 
functions extracted from the experimental facts in wide 
range of the variation of the scaled variable y. As seen 
from the mentioned Figures, the constructed normalized 
thermodynamic and transport functions show the scaling 
behavior over three decades in the normalized variable. 




B(T) 

FIG. 9: Temperature versus magnetic field T — B phase di- 
agram for YbRri2Si2. Solid circles represent the boundary 
between AF and NFL states. The solid squares denote the 
boundary of the NFL and LFL regime [Isl. |2& [3l| shown by 
the dotted line which is approximated by y/B — B c o [13 ]. Di- 
amonds mark the maximums T M of C/T m shown in Fig. 
[4j The dash-dot line is approximated by Tm oc a(B — B c o), a 
is a fitting parameter, see Eq. fl9l). Triangles along the solid 
line denote T in f in LMR [la . |20| ] sown in Fig. [8] the solid 
line represents the function T in f oc b(B ~ B c o), b is a fitting 
parameter, see Eq. (|10[) . 

Fig. IH1 reports T in f and Tm versus B depicted by the 
solid and dash-dotted lines, respectively. The bound- 
ary between the NFL and LFL regimes is shown by the 
dashed line, and AF marks the antiferromagnetic state. 
The corresponding data are taken from Ref. (l8l - [20l l3l| . 
It is seen that our calculations are in good agreement 
with the facts [l7|- In Fig. |H1 the solid and dash-dotted 
lines corresponding to the functions Ti„/ and Tm, respec- 
tively, represent the positions of the kinks separating the 
energy scales in C and M reported in Ref. [l8|, [2(J ■ It is 
seen that our calculations are in accord with facts, and 
we conclude that the energy scales are reproduced by 
Eqs. © and (fTU)) and related to the peculiar points of 
the normalized effective mass M^. The points are the in- 
flection point Ti n f and the maximum point Tm at which 
the transition region is located. These are shown by the 
arrows in the inset to Fig. [2j 

At B ->• B c0 both T inf -» and T M -> 0, thus the 
LFL and the transition regimes of both C/T and M 



V. UNIVERSAL BEHAVIOR OF 
TWO-DIMENSIONAL 3 He AT LOW 
TEMPERATURES 

The bulk liquid 3 He is historically the first object, 
to which the Landau Fermi-liquid (LFL) theory had 
been applied [26j . This substance, being an intrinsically 
isotropic Fermi-liquid with negligible spin-orbit interac- 
tion is ideal to test the LFL theory. It is remarkable 
that the same 3 He becomes the first 2D homogeneous 
Fermi-liquid in which the NFL behavior has been de- 
tected |39l - l4l| . 2D 3 He has a very important feature: a 
change of the number density x of 3 He film drives it to- 
wards QCP at which the quasiparticle effective mass M* 
diverges f39l - l4lj ]. This peculiarity permits to plot the 
experimental temperature-density phase diagram, which 
can be directly compared with the theoretical phase di- 
agram shown in Fig. [T] As a result, 2D 3 He becomes 
an ideal system to test a theory describing the NFL be- 
havior. Namely, the neutral atoms of 3 He are fermions 
interacting with each other by Van-der-Vaals forces with 
strong hardcore repulsion and a weakly attractive tail. 
The different character of inter-particle interaction along 
with the fact, that the mass of 3 He atom is 3 orders of 
magnitude larger than that of an electron, makes 3 He sys- 
tems to have drastically different properties than those of 
HF metals. Because of this difference nobody can be sure 
that the macroscopic physical properties of these systems 
will be more or less similar to each other at their QCP. 

In this Section we show that despite of very different 
microscopic nature of 2D 3 He and 3D HF metals, their 
main universal features at their QCP are the same, being 
dictated by the extended quasiparticles paradigm. Our 
analysis of the experimental measurements has shown 
that the behavior of 2D 3 He is quite similar to that of 
HF compounds with various ground state magnetic prop- 
erties. Namely, we demonstrate that the main universal 
features of 3 He experimental T-x phase diagram resemble 
those in HF metals and can be well captured utilizing the 
notion of FCQPT embracing the extended quasiparticles 
paradigm and thus deriving NFL properties of above sys- 
tems from the paradigm. We also show that the universal 
behavior of the effective mass of 2D 3 He coincides with 
that observed in HF metals. 
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A. The temperature-number density phase 
diagram of 2D 3 He 

As we seen in Section I, at QCP x = x c the effective 
mass diverges at T = and the leading term of this 
divergence given by Eq. ((4]) reads 



M*(x) 
M 



A 



B 



1 



x 

■ ' r: 



(14) 



Equation ([14]) is valid in both 3D and 2D cases, while 
the values of factors A and B depend on dimensionality 
and inter- particle interaction [l2| . At x > x c the fermion 
condensation takes place. Here we confine ourselves to 
the case x < x c . 

Equation (|14|) shows that the maximum value of the 
effective mass M\. ; oc 1/(1 — z) and it follows from ([6]) 
that Mm oc T~ 2 / 3 . As a result, we obtain that T M 
at which the effective mass reaches its maximum value 



Mlj oc T 2 / 3 is given by 



T M cc(l-z) 3 ^. 



(15) 



We note that obtained results are in agreement with nu- 
merical calculations 113, UM ■ 
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z=x/x c 

FIG. 10: The temperature-number density phase diagram 
of 2D 3 He. The part for z < 1 corresponds to HF be- 
havior divided into the LFL and NFL parts by the line 
Tm{z) oc (1 — z) s / 2 , where Tm is the effective mass max- 
imum temperature. The exponent 3/2 = 1.5 coming from 
Eq. (|15[) is in good agreement with the experimental value 
of 1.7 ± 0.1 The dependence M*[z) oc (1 - z) _1 shown 
by the dashed line points out QCP taking place at z = 1. 
The regime for z > 1 consists of the LFL piece (the shadowed 
region, beginning in the intervening phase z < 1 [3{J, which is 
due to the substrate inhomogeneities, see text) and the NFL 
regime at higher temperatures. 

In Fig. [TU1 we show the phase diagram of 2D 3 He in the 
variables T - z (see Eq. (fT4]0 . For the sake of comparison 
the plot of the effective mass versus z is shown by dashed 
line. The dependence M*(z) oc (1 — z)^ 1 demonstrates 
that the effective mass diverges at QCP with z = 1 in 
accordance with the general phase diagram displayed in 
Fig. [T] The part of the diagram where z < 1 corresponds 
to HF behavior and consists of LFL and NFL parts, di- 
vided by the line Tm{z) oc (1 — z) 3 ! 2 . We draw attention 



here, that our exponents 1 (see Eq. (|T4"|) ) and 3/2 = 1.5 
(see Eq. (fTS"|) ) are in good agreement with these from 
Ref. [39| ■ The good coincidence between the theoretical 
and experimental exponents speaks in favor of realization 
of our FCQPT scenario in the NFL behavior of both 2D 
3 He and HF metals as former system is in great detail 
similar to them. 

The regime for z > 1 located above the quantum crit- 
ical line, see Figs. HTHandfTl consists of low-temperature 
LFL piece, (shown in Fig. [TU] by shadowed region, be- 
ginning in the intervening phase z < 1 [39]) and NFL 
regime at higher temperatures. The former LFL piece 
is related to the peculiarities of substrate on which 2D 
3 He film is placed. Namely, it is related to weak sub- 
strate heterogeneity (steps and edges on its surface) so 
that quasiparticles, being localized (pinned) on it, give 
rise to the LFL behavior [3^, H(| ■ That is the peculiari- 
ties of the substrate eliminate the degeneracy generated 
by the FC state taking place at z > 1 in the same way 
as the AF state does in the case of YbRh 2 Si2, see Fig. 
At elevated temperatures, the competition between ther- 
mal and pinning energies returns the system back to the 
unpinned state restoring the NFL behavior. As HF met- 
als do not have a substrate, the LFL behavior is induced 
by the AF state lifting the degeneracy. At elevated tem- 
peratures, this state is destroyed and exhibits the NFL 
behavior, as it is shown in Fig. [5] If the AF state were 
absent and some disorder (like point defects, dislocations 
etc) were present in the lattice a rather thin LFL piece 
could take place at low temperatures. 



B. NFL behavior of 2D 3 He versus that of HF 
metals 



M/M 



M/M 




FIG. 11: The dependence of the effective mass M* (z) on the 
dimensionless density 1 — 2 = 1 — x/x c . Experimental data 
from Ref. 14111 are shown by circles and squares and those 
from Ref. [3y] are shown by triangles. The effective mass is 
fitted as M*(z)/M otA + B/(l - z) (see Eq. fli)). while the 
reciprocal one as M/M*(z) oc Aiz, where A, B and Ai are 
constants. 

As we have seen above, M*(T) can be measured in 
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experiments on strongly correlated Fermi-systems. For 
example, M*(T) cx C(T)/T cx S(T)/T cx M (T) cx X (T) 
where C{T) is the specific heat, S(T) — entropy, M (T) 
— magnetization and x(T) — AC magnetic susceptibil- 
ity. If the measurements are performed at fixed x then, 
as it follows from Eq. ([Jj), the effective mass reaches the 
maximum at T = Tm- Upon normalizing both M*(T) 
by its peak value at each x and the temperature by Tm, 
we see from Eq. ([7]) that in the case of 2D 3 He all the 
curves also merge a into single one, demonstrating a scal- 
ing behavior. 




Normalized temperature 

FIG. 12: The normalized effective mass M"^ as a function 
of the normalized temperature T '/Tm at densities shown in 
the left bottom corner. The behavior of M% is extracted 
from experimental data obtained in 2D 3 He [4Q| and 3D HF 
compounds with different magnetic ground states such as 
CeRu2Si2 and CePdi- x Rh x [U 5^, fitted by the solid curve 
given by (UJ. 

In Fig. II H we report the experimental values of ef- 
fective mass M*(z) obtained by the measurements on 
3 He monolayer [41| . These measurements, in coinci- 
dence with those from Ref. [3!|, show the divergence 
of the effective mass at x = x c . To show that our FC- 
QPT approach is able to describe the above data, we 
represent the fit of M*(z) by the fractional expression 
M*(z)/M cx A + -8/(1 — z) and the reciprocal effective 
mass by the linear fit M/M*(z) cx A\Z. We note here, 
that the linear fit has been used to describe the exper- 
imental data for a bilayer of 3 He [39| and we use this 
function here for the sake of illustration. It is seen from 
Fig. [TT] that the data [3!| ( 3 He bilayer) can be equally 
well approximated by both linear and fractional func- 
tions, while the data [4l| cannot. For instance, both 
fitting functions give for the critical density in bilayer 
x c R3 9.8 ran -2 , while for monolayer [4l| these values are 
different: x c = 5.56 for a linear fit and x c = 5.15 for a 
fractional fit. It is seen from Fig. [11] that a linear fit is 
unable to properly describe the experiment (4lj at small 
1 — z (i.e. near x = x c ), while the fractional fit describes 
the experiment very well. This means that more detailed 
measurements are necessary in the vicinity x = x c . 

We now apply the universal dependence ([7]) to fit the 
experiment not only in 2D 3 He but also in 3D HF metals. 
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FIG. 13: The dependence M^(T/Tm) at densities shown in 
the left down corner. The behavior Afjy is extracted from 
experimental data for C(T)/T in 2D 3 He [HI and for the 
magnetization Mo in 2D 3 He [3^]. The solid curve shows the 
universal function, see the caption to Fig. [T2l 

Mff(y) extracted from the entropy S(T)/T and magne- 
tization Mo measurements on the 3 He film (4(| at dif- 
ferent densities x < x c is reported in Fig. [12] In the 
same figure, the data extracted from the heat capacity 
of ferromagnet CePd0.2R.h0. 8 [HJ and the AC magnetic 
susceptibility of paramagnet CeRu2Si2 [HJ are plotted 
for different magnetic fields. It is seen that the universal 
behavior of the normalized effective mass given by Eq. 
([7]) and shown by the solid curve is in accord with the 
experimental facts. All 2D 3 He substances are located at 
x < x c , where the system progressively disrupts its LFL 
behavior at elevated temperatures. In that case the con- 
trol parameter, driving the system towards its QCP x c is 
merely the number density x. It is seen that the behavior 
oiM* N {y), extracted from S(T)/T and M of 2D 3 He (the 
entropy S(T) is reported in Fig. S8 A of Ref. 40]) looks 
very much like that of 3D HF compounds. As we shall 
see from Fig. [T4l below, the amplitude and positions of 
the maxima of magnetization Mq(T) and S(T)/T in 2D 
3 He follow nicely Eqs. (TT4")) and (IT5"1) . We conclude that 
Eq. allows for the reduction of a 4D function describ- 
ing the effective mass to a function of a single variable. 
Indeed, the effective mass depends on the magnetic field, 
temperature, number density and composition so that all 
these parameters can be merged in the single variable by 
means of interpolating function like Eq. ([Jj. 

The attempt to fit the available experimental data for 
C(T)/T in 2D 3 He by the universal function M* N (y) 
is reported in Fig. [JjJ] Here, the data extracted from heat 
capacity C(T)/T for the 3 He monolayer [4l[ and magne- 
tization Mo for the bilayer [3^], are reported. It is seen 
that the effective mass extracted from these thermody- 
namic quantities can be well described by the universal 
interpolation formula ([Jj. We note the qualitative and 
quantitative similarity between the double layer [3^ | and 
the monolayer [4l[ of 3 He as seen from Fig. [JjJJ 

In the left panel of Fig. [HI we show the density de- 
pendence of Tm, extracted from measurements of the 
magnetization Mq(T) of the 3 He bilayer [39]. The peak 
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FIG. 14: Left panel, the peak temperatures Tm and the peak 
values M m ax extracted from measurements of the magnetiza- 
tion Mq in 3 He 39]. Right panel shows Tm and the peak 
values (5/T)max extracted from measurements of S(T)/T in 
3 He We approximate Tm oc (1 - z) 3/2 and (5 , /T) max oc 
l^aA/fl-z). 

temperature is fitted by Eq. (fTS"|). In the same Figure, 
we have also reported the maximal magnetization Af max . 
It is seen that M max is well described by the expression 
M max oc (5/T) max oc (1 - z)~\ see Eq.' (HJ). The right 
panel of Fig. [14] reports the peak temperature Tm and 
the maximal entropy (<S'/T) max versus the number den- 
sity x. They are extracted from the measurements of 
S{T)/T on the 3 He bilayer 0. The fact that both the 
left and right panels exhibit the same behavior of the 
curves shows once more that there are indeed the quasi- 
particles, determining the thermodynamic behavior of 2D 
3 He near its QCP related to FCQPT. 

VI. SUMMARY 

We have analyzed the non- Fermi liquid behavior of the 
heavy fermion metals, and showed that extended quasi- 



particles paradigm is strongly valid, while the depen- 
dence of the effective mass on temperature, number den- 
sity and applied magnetic fields gives rise to the NFL be- 
havior. We have demonstrated that our theoretical study 
of the heat capacity, magnetization, longitudinal mag- 
netoresistance and magnetic entropy are in good agree- 
ment with the outstanding recent facts collected on the 
HF metal YbRli2Si2. Our normalization procedure has 
allowed us to construct the scaled thermodynamic and 
transport properties in wide range of the variation of the 
scaled variable. For YbRh 2 Si2 the constructed thermo- 
dynamic and transport functions show the scaling behav- 
ior over three decades in the normalized variable. The 
energy scales in these functions are also explained. 

We have described the diverse experimental facts re- 
lated to temperature and number density dependencies of 
the thermodynamic characteristics of 2D 3 He by a single 
universal function of one argument. The above univer- 
sal behavior is also inherent to HF metals with different 
magnetic ground states. We obtain the marvelous coin- 
cidence with experiment in the framework of our theory. 
Moreover, these data could be obtained for 2D 3 He only 
and thus they were inaccessible for analysis in HF met- 
als. This fact also shows the universality of our approach. 
Thus we have shown that bringing the experimental data 
collected on different strongly correlated Fermi-systems 
to the above form immediately reveals their universal 
scaling behavior. Thus, the theory of fermion condensa- 
tion quantum phase transition, preserving the extended 
quasiparticles paradigm and intimately related to unlim- 
ited growth of the effective mass as a function of tem- 
perature, magnetic field etc, is capable of describing the 
strongly correlated Fermi systems. 
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